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Unfortunately,  previous  research  did  not  show  generally  how  to  determine  this  parameter.  We  present 
general  theoretical  considerationt  for  its  numerical  noniterative  determination.  This  general  theory  it 
igplied  to  an  important  clast  of  attrition.rate  coefficients  (offset  powerMate  coefficients).  Our  results  allow 
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James  G.  TAYumt  and  Gerald  G.  Brown^ 

Depirtmenu  of  Operations  Research  and  Computer  Science,  Naval  Postfiaduate  School,  Monterey, 

California 

Seapc  aad  purpam— This  paper  presents  new  computational  methods  that  facilitate  digital<omputer 
analysis  of  some  imp.  tant  military  operations  research  problems.  Lanchester-type  combat  models(l]  are 
deterministic  differential.equation  mo^ls  of  combat  attritioo  in  which  the  state  variables  are  typically  the 
numbers  of  the  different  weapon-system  types.  Even  though  combat  between  two  military  forces  is  a 
complex  random  process,  such  deterministic  combat  models  are  commonly  used  for  computational  reasons 
in  ddense-planainc  studies,  for  example,  to  assess  the  relative  importance  of  various  weapon-system  and 
force-level  parameters.  A so-caHed  attrition-rate  coeflcient  in  such  a combat  model  represents  the  fire 
elfectivenest  of  a particular  weapon-system  type  against  a particular  target  type,  and  time-dependent 
attrition-rate  coefficientt  are  used  to  model  temp^  variatkms  in  fire  effectiveness  when,  for  example,  the 
range  between  firers  and  targeu  changes  appreciably  dining  battle.  For  such  a variable-coeflicient 
Lanchester-type  combat  model  that  is  a generalization  of  Lancbester’t(21  clastic  “square-law”  model,  we 
present  a simple  numerical  procedure  for  determining  the  so-called  parity-conditioa  parameter,  which  it 
“the  enemy  force  equivalent  of  a friendly  force  of  unit  strength”  and  may  be  used  to  predict  battle  outcome 
in  specific  engagememt.  These  resuhs  allow  one  not  only  to  predict  battle  outcome  but  also  to  tradeoff 
quality  vs  quantity  of  two  opposing  weapon  systems. 

AhMract— This  paper  presents  a simple  numerical  procedure  for  determining  the  parity-condition  parameter 
for  determinittic  Lanchester-type  combat  between  two  homogeneous  forces.  Deterministic  differential- 
equation  combat  modelt  ate  commonly  used  in  parametric  studiet  for  computational  reasons,  since  they 
give  essentially  the  same  resuhs  for  tk  mean  course  of  combat  as  do  corresponding  stochastic  attrition 
models.  The  combat  studied  in  this  paper  it  modelled  by  Lanchester-type  equations  of  modern  warfare  with 
time-dependent  attrition-rate  coefficients.  Previous  research  hat  generalized  Lanchester's  clastic  “square 
law"  to  such  variable-coefBcient  combat.  It  has  shown  that  the  prediction  of  battle  outcome  (in  particular, 
force  annihilation)  without  having  to  spend  the  time  and  effort  of  computing  force-level  trajectories 
depends  on  a single  parameter,  the  so-called  parity-condition  parameter,  which  it  “the  enemy  force 
equivalent  of  a friendly  force  of  unit  strength”  and  depends  on  only  the  attritioo-rate  coefficients. 
Unfortunately,  previous  research  did  not  show  generally  how  to  determine  this  parameter.  We  present 
general  theoretical  consideratioos  for  its  numerical  noniterative  determination.  This  general  theory  is 
applied  to  an  important  class  of  attrition-rale  coefficients  (offset  power-rate  coefficients).  Oir  results  allow 
one  to  study  such  variable-coefficient  combat  models  almost  as  easily  and  tbortM^y  as  Lanchester's 
classic  constant-coefficient  model. 


I.  INTRODUCTION 

As  a consequence  of  pioneering  work  by  F.  W.  Lanchester  done  about  the  time  of  World  War  I, 
military  analysts  have  used  simplified  deterministic  differential-equation  models  to  develop 
insighu  into  the  dynamics  of  combat  from  about  the  end  of  World  War  II  (see,  for  example, 

[1, 3-11].  Such  deterministic  models  have  been  widely  used  because  (among  other  reasons)  the 

*Tliit  research  was  partiaOy  supported  by  the  Office  'of  Naval  Research  and  by  the  U.S.  Aimy  Research  Office, 

IXiiliaffl,  North  Carolina,  under  RAD  Project  No.  1LI6I102H57-05  Math. _» 

tJames  G.  Taylor  is  an  Associate  Professor  of  uperatiou  Research  at  the  ulsTNaval  Postgraduate  School  (NPS).  He 
received  B.S.,  M.S.,  and  Ph.D.  degrees  in  Petroleum  Engineering  from  Stanford  University.  He  is  the  author  of  over 
twenty  puMicatioos  on  various  aspeett  of  operatioos  research  aad  was  awarded  (jointly  with  S.  Parry)  the  1975  MAS  Prize 
by  the  Military  Applications  Section  (MAS)  of  the  Operalioiis  Research  Society  of  America  (ORSA).  He  is  the  co-author 
with  Professor  Gerald  G.  Brown  of  severM  papers  on  Lanchester-type  models  aad  optimal  control  theory.  Professor 
Taylor's  recent  research  hu  been  on  the  Lanchester  theory  of  combat  and  generalized  control  theory.  He  serves  u an 
Associate  Editor  of  both  Opemtkmt  Research  and  Naval  Kattarck  LogMet  Qaantriy  and  also  u a Councilman  of  the 
Military  Applications  Section  of  ORSA. 
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corresponding  stochastic  formulations  are  for  all  practical  purposes  analytically  intractable  (see 
Note  1 on  p.  65  of  [12]).  The  advent  of  the  modern  high-speed  digital  computer  has  made  feasible 
the  development  and  use  of  quite  complicated  versions  of  such  Lanchester-type*  models  as 
practical  defense  planning  tools[13].  Thus,  today  militarily  realistic  computer-based  Lanchester- 
type  models  of  quite  complex  combat  systems  have  been  developed.  Such  models  currently  exist 
for  almost  the  entire  spectrum  of  combat  operations,  from  combat  between  battalion-sized  ( 14]  and 
division-sized  [IS]  units  to  theater-level  operations  [16, 17]. 

A simple  combat  model,  however,  may  yield  a clearer  understanding  of  important  inter- 
relationships that  are  difficult  to  perceive  in  a more  complex  model,  and  such  insists  can 
provide  valuable  guidance  for  more  detailed  computerized  investigations  (see  [4, 1 1]).  In  this 
paper  we  present  a new  important  numerical  procedure  that  facilitates  parametric  analysis  (in 
particular,  the  parametric  examination  of  force-annihilation  prediction)  of  battle  outcomes  for 
such  simplified  Lanchester-type  models  of  combat  between  two  homogeneous  forces  with 
temporal  variations  in  each  side’s  fire  effectiveness.  Previously,  such  battle-outcome  in- 
formation couk)  only  be  readily  obtained  from  constant-coefficient  models,  and  S. 
Bonder  [3, 18, 19]  has  emphasized  the  deficiencies  of  constant-coefficient  models  (see  Section  3 
below).  These  results  are  not  only  significant  in  their  own  right  but  are  also  useful  in  the 
quantitative  analysis  of  time-sequential  combat  strategies  (see,  for  example,  [20, 21]). 

It  is  important  for  the  military  operations  analyst  to  have  a clear  understanding  of  how 
force-level  and  weapon-system-performance  parameters  interact  to  determine  a battle’s 
outcome.  Such  knowledge  is  particularly  useful  in  weapon-system  and  force-level  planning 
activities  for  defense  planning  (especially  since  one  frequently  uses  models  that  are  so 
complicated  that  trends  are  not  directly  discernible  without  extensive  (and  costly)  computer 
runs).  S.  Bonder’s[3,4, 19]  pioneering  work  on  methodology  for  the  evaluation  of  military 
systems  (particularly  mobile  systems  such  as  tanks,  mechanized  infantry  combat  vehicles,  etc.) 
provides  a motivation  for  interest  in  variable-coefficient,  deterministic,  Lanchester-type  combat 
models  such  as  we  consider  in  this  paper.  He  has  stressed  (see  pp.  30-31  of  [4])  the  importance 
of  analytical  solutions  to  such  models  for  developing  insights  into  the  dynamics  of  combat  by 
portraying  the  relation  between  various  factors  in  the  combat  attrition  process  and  the  surviving 
numbers  of  forces  and  for  facilitating  sensitivity  and  other  parametric  analyses  (see  [22]). 
Unfortunately,  as  work  by  Bonder  and  Farrell[4]  and  Taylor[12,23]  shows,  the  analytical  (i.e. 
infinite  series)  solution  to  variable-coefficient  equations  generally  by  itself  (i.e.  without  expli- 
citly computing  force-level  trajectories)  provides  little  information  about  battle  outcome 
because  of  its  complexity.  Therefore,  one  must  seek  new  ways  for  developing  insights. 

Taylor  and  Comstock[7]  have  given  results  that  allow  one  to  predict  battle  outcome  (in 
particular,  force  annihilationt)  in  theory  without  having  to  spend  the  time  and  effort  of 
computing  force-level  trajectories.  To  be  computationally  practical,  however,  their  results 
require  the  determination  of  the  so-called  parity-condition  parameter  (“the  enemy  force 
equivalent  of  a friendly  force  of  unit  strength"),  which  depends  on  only  the  model’s  attrition- 
rate  coefficients.  They  analytically  determine  the  parity-condition  parameter  for  power  attrition- 
rate  coefficients  with  “no  offset”,  which  allow  one  to  model  combat  between  two  weapon 
systems  with  the  same  minimum  effective  range  but  different  range  dependencies  for  each 
system’s  fire  effectiveness  (see  also  [6]).  It  is  the  purpose  of  this  paper  to  show  how  to 
determine  the  parity-condition  parameter  in  other  cases,  in  particular  for  power  attrition-rate 
coefficients  with  “positive  offset”,  which  allow  one  to  model  such  combat  between  weapon 
systems  with  different  maximum  effective  ranges.  Our  results  allow  one  to  study  in  general  such 
variable-coefficient  combat  models  almost  as  easily  and  thoroughly  as  Lanchester’s  classic 
constant-coefficient  model. 

The  organization  of  this  paper  is  as  foUows.  We  first  review  Lanchester-type  equations  of 
modem  warfare,  especially  variable-coefficient  formulations.  Next  we  review  force-anni- 
hilation-prediction  conditions  for  such  models  and  show  how  to  use  knowledge  about  the 

*Mio  frequently  ceBed  ddferetiliti  modeli  of  co«httfl5]. 

tBonder  tad  Honi((5]  point  o«,  however,  that  force  innihilttion  aiay  not  ahrtyt  be  the  beet  aiMrioa  for  iTriiiatim 
mflitiry  opcraiioat.  See  pp.  192-242  of  Bonder  ind  PatreDldl  tor  t dctoaed  Laacbouer-typo  anotyih  of  an  atUck  tcanorie 
for  whicb  other  “end  of  battle  conditioni'’  phy  the  principal  role.  Neverlboleaa.  ilia  of  conildwibli  lalweat  (im»rioly  ter 
developins  intishu  into  the  dyntmici  of  combat)  to  be  able  to  early  predict  the  ocewronco  of  tareo  aaaMtlien. 
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parity-condition  parameter  for  one  set  of  attrition-rate  coefficients  to  numerically  determine  it 
in  related  cases  of  interest.  This  general  theory  is  then  applied  to  the  important  case  of  offset 
power  attrition-rate  coefficients,  with  detailed  numerical  examples  being  given. 

2.  LANCHESTER'S  CLASSIC  COMBAT  FORMULATION 
F.  W.  Lanchester[2]  (see  also  p.  45  of  [12])  hypothesized  in  1914  that  combat  between  two 
military  forces  could  be  modelled  by 

dxldt  = -ay,  dy/d(>-hx,  (1) 

with  initial  conditions 

x(t*0)=xo,  y(t»0)«y^  (2)  I 

where  t = 0 denotes  the  time  at  which  the  battle  begins,  x(()  and  y(t)  denote  the  numbers  of  X ^ 

and  Y at  time  t,  and  (for  a particular  battle)  a and  b are  nonnegative  constants  which  are  today 
called  Lanchester  attrition-rate  coefficienU  and  represent  each  side’s  fire  effectiveness.  The 
equations  (1)  are  only  valid  for  x,  y > 0.  For  example,  the  first  becomes  dx/dt  = 0 when  x = 0. 

Because  of  Lanchester’s  pioneering  work  [2],  we  will  refer  to  any  differential-equation 
model  of  combat  attrition  as  a Lanchrster-type  combat  model  or  as  a system  of  Lanchester- 
type  differential  equations  (or  sometimes  simply  as  Lanchester-type  equations).  In  particular, 
we  will  refer  to  (1)  as  (constant-coeiicient)  LoHckesler-type  equations  of  modem  warfare. 

Other  forms  of  Lanchester-type  equaiioas  appear  in  the  literature [1,23],  but  we  will  not 
consider  these  here.  Various  sets  of  physical  circumstances  have  been  hypothesized  to  yield 
them:  for  example:  (a)  both  sides  use  aimed  fire  and  target  acquisition  times  are  constant[10]; 
or  (b)  both  sides  use  area  fire  and  a constant  density  defense  (see  p.  345  of  [1]). 

From  (1)  Lanchester  deduced  his  famous  square  law 

bix^-xHt))^a{yo^-y\t)).  (3)  | 

Consider  now  a battle  terminated*  by  either  force  level  reaching  a given  “breakpoint”:  for 
example,  Y wins  when  X/  » x(tf)  - Xbp  * /x*'x.  but  yf>ytp~  /r^'y©.  where  //,  x/,  y/  denote 
final  values  and  Xgp  denotes  X's  breakpoint  which  is  a given  fraction  /x*'  of  his  ini^  strength. 

It  follows  from  (3)  that 


Y will  win  if  and  only  if  ^ (mi-(/x"^)’ 

which  for  a fight-to-the-finish  (i.e.  /x*'’  = /v*^  = 0)  becomes  the  classic  result 

Y will  win  a fight-to-the-finish  if  an  only  ify<-^  (5) 

Unfortunately,  no  relationship  similar  to  (3)  holds  in  general  for  variable  attrition-rate 
coefficients  except  when  a(t)lb{t)  - constant  (see  p.  48  of  [12]).  This  paper,  nevertheless,  shows 
how  (5)  generalizes  in  these  cases,  but  so  far  we  iuve  not  been  able  to  generalize  (4).  Recalling 
that  the  time  history  of  the  X force  level  is  given  by 

x(l)  - xo  cosh  \/(ab)t  “ »inli  y/(ob)t,  (6) 

we  see  that  the  battle  trajectories  depend  on  the  two  weapon-system-performance  parameters: 
(I)  the  intensity  of  combat  \/ab,  and  (11)  the  relative  fire  effectiveness  alb.  Only  the  relative  fire 
effectiveness,  however,  determines  the  battle’s  outcome  [see  (4)  and  (5)  above]. 

The  ■odeUai  of  bank  tcfariBedoa  it  ■ proMw  ena  is  ceasiwpenfy  Meate  phaaim  Madiee  (lee  pa.  S24-S2S  of 
HI). 
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3.  VARIABLE  ATTRITION-RATE  COEFFICIENTS 
Bonder[I81  has  pointed  out  that  in  many  cases  (for  example,  in  the  case  of  mobile  weapon 
systems)  the  validity  of  the  assumption  of  constant  attrition-rate  coefficients  is  open  to  question 
(see  also  [3, 4, 19]).  Thus,  we  consider 


dxidt  - - a(f)y,  dy/dr  = - h(/)x,  (7) 

where  a(t)  and  b(.t)  denote  time-dependent  attrition-rate  coefficients.  We  assume  that  aU)  and 
b(l)  are  defined,  positive,  and  continuous  for  to<t<  + ’*>  with  to^O.  For  convenience,  we 
introduce  the  notation  that  aU)BL(t^T)  means  exists  (and  is  given  by  a finite 

quantity).  From  our  assumptions  about  a(t)  and  b(0,  it  follows  that  a(()€  L((o.  T)  implies  that 
/,[  a(t)  dt  = +oo.  We  also  assume  that  a(r),  b(t)  e Wo,  T)  for  any  finite  T.  We  further  take  a(t) 
and  b(t)  to  be  given  in  the  form  a(t) » b(t)  M(0,  where  k,,  ko  are  positive  constants 
chosen  so  that  aU)lbU)  = kjko  when  g(t)  > h(().  Analogous  to  the  constant-coefficient  case 
[see  discussion  after  (6)],  we  have  the  two  weapon-system-performance  parameters:  (I)  the 
intensity  of  combat,  I(t)  = \/{a(t)b(t)y,  and  (II)  the  relative  fire  effectiveness,  R(t)  > a(()/h(t). 
We  accordingly  introduce  the  combat-intensity  parameter  A/  and  the  relative-fire-ellectiveness 
parameter  Ar  defined  by 

A/  = V(M»)  and  A*  = kjk».  (8) 

Two  significant  developments  in  the  Lanchester  theory  of  combat  during  the  1960s  were  the 
development  of  methodology  for  (a)  the  prediction  of  Lanchester  attrition-rate  coefficients  from 
weapon-system-performance  data  by  S.  Bonder[18,24]  and  (b)  the  (maximum  likelihood) 
estimation  of  such  coefficients  from  Monte  Carlo  simulation  output  by  G.  Clark[2S].  Both  these 
developments  and  others  (see  [12]  for  further  references)  have  generated  interest  in  the  model 
(7)  and  facilitated  its  application  (and  that  of  its  generalization  to  combat  between  hetero- 
geneous forces  [4)  to  defense  planning  studies. 

A large  class  of  tactical  situations  of  interest  can  be  modelled  with  the  following  general 
power  attrition-rate  coefficients  [4, 7, 12] 


k.(f  + Cr,  and 


= Mt  + C + A)', 


where  A,  C 2 0.  We  will  call  A the  offset  parameter,  since  it  allows  us  to  model  (with  ii,v^0) 
battles  between  weapon  systems  with  different  maximum  effective  ranges.  We  will  call  C the 
starting  parameter,  since  it  allows  us  to  model  (again  with  ^ 0)  battles  that  b^in  within  the 
minimum  of  the  maximum  effective  ranges  of  the  two  systems.  For  example,  let  us  consider 
Bonder's  [3, 19]  model  of  a constant-speed  attack  on  a static  defensive  position  (see  also 
[12, 23]).  Then  we  have 

dx/dl  * - o(r)y,  dy/dr  = - ^(r)x,  (10) 

where  r(0  ■ ff«-  M denotes  the  distance  (range)  between  the  two  opposing  forces,  Ro  denotes 
the  bank's  opcnmg  range,  v > 0 denotes  the  constant  attack  speed. 


/»  fO  torr^Rt, 

*^^^“la,(l-f/l«.r  forOsrsJL 

is  20,  and  X,  denotes  the  maximum  effective  range  of  Y't  weapon  system.  Similarly  for  fi(r), 
with  exponent  r 2 0.  In  (11)  the  parameter  /i  allows  us  to  model  the  range  dependence  of  Y’s 
fire  effectiveness  (see  Rg.  1).  The  offset  and  starting  parameters  are  given  by 

A^iRp-RoVv,  and  C-(X.-XJ/o.  (12) 

and  the  assumption  A,C20  nplks  that  2 Jf.  2 Jl».  Ftom  considering  (12)  and  Fig.  2,  the 
reader  should  have  no  troubk  understanding  our  terminology  for  A and  C. 


.i  liUi  V. 


I 

)■ 


Lanchester-type  equatioiu 


Attrition- 

Rot* 

Cooffieiont 
a (A)  0.4 


0.2  + 


iOOO  ■ iSOO 
Rong*  /I  (motors) 

Fig.  I.  Dependence  the  aurition-rate  coefficient  a(r)  on  tlie  exponent  ft  with  maximum  effective  range  c 
the  weapon  system  and  kill  capability  at  zero  range  held  constant.  (Notes:  (I)  the  maximum  effective  ranged 
of  the  system  is  denoted  as  X.  > 2000  m;  (2)  a(r  > 0)  • oo  > 0.6X  casualties/(unit  time  x number  of  Y units) 
denotes  the  F-force  weapon-system  kill  rate  at  zero  force  separation  (denoted  here  as  range);  (3)  the 
opening  range  of  battle  is  denoted  as  Xq  > 1250  m and  (as  shown)  Xo  < X..) 


f 
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I 
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Fig.  2.  Explanation  of  offset  parameter  A and  star^  parameter  C for  power  attrition-rate  coefficients 
■ndiling  constaM-specd  attack.  (Notes:  (1)  the  maximum  effective  ranges  of  the  two  weapon  systems  are 
denotnd  at  Xi  and  (2)  the  opening  range  of  battle  (i.c.  initial  separation  between  forces)  it  denoted  at 
Xs  and,  at  shosvn,  Xa  < Minimum  (Xn  Xf);  (3)  the  offset  parameter  n given  by  A - (J^  - X.)/e;  (4)  the 
starting  parameter  is  given  by  C - (X.  - Xo)/v.l 

The  time  history  of  the  A’  force  level,  i.e.  the  solution  x(f)  to  (7),  is  given  by  [12] 

x(t)  - Xo(Cy(0)Cx(f)  - Sr(0)Sx(/))  - yoV(A«MCx(0)5x(f ) - Sx(0)Cx(/)}.  (13) 

where  the  hyperbolic-like  generul  Lanchester  functions  (GLF)  Cx(f)  and  5x(f)  are  linearly 
independent  solutions  to  the  A force-level  equation 

with  initial  conditions 

Cx(t-t*)-l.  Sx(f-t*)-0. 

{[l/a(f)JdCx/df(f)},.„-0.  {(I/«(f)]d5x/df(f)},.^-I/V(A*). 
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where  to  denotes  the  la.'gest  finite  time  at  which  d(r)  or  b(t)  ceases  to  be  defined,  positive,  or 
continuous.  More  precisely,  to  ~ sup  {(||either  a(t)  or  b(t)  is  not  defined,  positive,  or  continuous 
for  some  finite  r 2 tj  = inf  {t||both  a(t)  and  b(t)  are  defined,  positive,  and  continuous  for  all 
finite  t > to}.  For  example,  to  ^ ~ C for  the  general  power  attrition-rate  coeflBcients  (9).  The  time 
history  of  the  F force  level  may  be  similarly  obtabed,  with  Cr(t)  and  Sr(0  being  analogously 
defined  for  the  corresponding  Y force-level  equation. 

For  the  numerical  determination  of  the  parity-condition  parameter,  it  is  convenient  to 
introduce  a new  independent  variable  s defined  by 


s^KXi  r g(a)A<r, 


where  the  parameter  K is  to  be  chosen  to  simplify  the  form  of  J(s)  given  by  (18).  We  denote 
5(r  ==  0)  as  5o,  and  then  so^O  if  and  only  if  to^O.  The  substitution  (16)  transforms  (14)  into  the 
normal  form 

^-/(j)x  = 0.  (17) 


. 1 

mw)}' 


and  t = t(s).  We  also  define  the  normal-form  hyperbolic-like  GLF  cx(s)  and  sx(s),  which 
satisfy  (17)  and  the  initial  conditions 


cx(*  = 0)»l. 


Jx(«=0)  = 0, 


dcx/d5(s  = 0)  = 0,  d5x/ds(s  = 0)  = 1. 


It  follows  that 


Cx(s(/))=Cx(f).  and  sxis(t))  = KSxU).  (20) 

4.  FORCE-ANNIHILATION-PREDICTION  CONDITIONS 
Recently,  Taylor  and  Comstock  [7]  have  generalized  the  constant-coeflScient  force-anni- 
hilation-prediction condition  (3)  to  battles  modelled  with  Lanchester-type  equations  of  modern 
warfare  with  time-dependent  attrition-rate  coefficients  (7).  In  some  sense  their  results  generalize 
Lanchester’s  famous  square  law  to  variable-coefficient  combat.  Taylor  and  Comstock  have 
shown  that  for  the  model  (7)  force-annihilation  prediction  involves  (besides  the  initial  force 
ratio  xjyo  and  a relative-fire-effectiveness  parameter  A*)  a single  parameter,  which  they  call  the 
parily-condUioH  parameter,  denoted  here  as  (?*.  Their  main  theoretical  result  is  stated  here  for 
the  reader’s  convenience  as  Theorem  1. 

Theorem  1 

Taylor  and  ComstockfT],  assume  that  either  tt(t)tEU0,+<*>)  or  h(()€L(0, -«-<>>).  Then  the 
X force  will  be  annihilated  in  finite  time  1/  and  only  if 

where  the  parity-condUion  parameter  Q*  is  unUiue  and  given  by 

ii«  m JL  (22) 

•"cx(o  <?•• 


Linchester-type  equations 


233 


I 

i 


Remark  1.  The  condition  (21)  means  that  neither  force  is  annihilated  (in  other  words,  the 
forces  are  of  “equal  fighting  strength”)  if  and  only  if 


yo 


Cx(0)-<?*5x(0) 

Q*Cy(0)-Sy(0) 


}■ 


(23) 


which  when  to  = 0,  simplifies  to 


Xq  V(Ak) 

yo  Q*  • 


(24) 


In  other  words,  the  above  equation  (23)  is  a condition  under  which  two  forces  are  “at  parity” 
with  one  another  (hence,  the  term  parity-condition  parameter).  Observing  the  special  case  (24), 
we  may  consider  the  parity-condition  parameter  to  be  “the  enemy  force  equivalent  of  a friendly 
of  unit  strength”  (see  also  equation  (17)  of  [7]). 

Remark  2.  We  also  have  lim,.^ {SyfO/Cj-lt)}  = (?*. 

Remark  3.  When  to-0,  (21)  simplifies:  X will  be  annihilated  in  finite  time  if  and  only  if 
JCo/yo<V(A«)/Q*. 

Remark  4.  The  result  (22)  suggests  a numerical  procedure  for  approximately  determining 
the  parity-condition  parameter  Q*:  we  may  approximate  the  parity-condition  parameter  Q*  by 
» ll{Sx(t)ICxU)h  where  t is  a “suitably  la^”  value  of  t.  In  other  words,  we  may  estimate 
Q*  simply  by  picking  a large  value  for  t (we  denote  this  selected  large  value  by  (),  computing 
Sx(t)  ^ Cx(0.  and  then  forming  their  ratio.  Our  estimate  for  Q*  is  tlwn  given  by  (^  = 
l/{Sx(t)/Cx(t)}.  The  only  problem  is  that  we  do  not  know  how  large  to  take  t for  “satisfactory” 
estimation  of  Q*:  There  is  an  estimation  error,  £(r)  > Q*  - (^f),  which  depends  monotonically 
on  t,  and  a priori  we  do  not  know  how  large  this  error  is.  The  present  paper  develops  a bound 
on  the  magnitude  of  this  error,  and  our  new  error  estimate  allows  the  goodness  of  ap- 
proximation to  be  easily  evaluated  in  many  cases  of  interest. 

We  may  also  determine  the  parity-condition  parameter  with  the  normal-form  hyperbolic-like  ’ 
GLF,  since  lim.^{sx(s)/cx(s)}=  I/Z*>  where  Z*  is  called  the  modified  parity- 
condition  parameter.  In  fact,  we  wfll  find  it  more  convenient  to  do  so.  With  this  in  mind,  let  us 
introduce  the  T-functions  Cv(s)  and  sy(s)  (corresponding  to  Cx(s)  and  Sx(s)]  defined  by 


dcy/ds  > /(5)sx,  dsy/ds  “ y(s)Cx,  (25) 

with  initial  conditions 

Cy(s -0)-l,  Sy(s*0)»0. 

It  follows  that  Cy(s)  and  sr(s)  are  linearly  independent  solutions  to  the  modified  Y equation 


ds{/(s)ds} 


(26) 


and 


CrUU))  - CrU).  ir(s(l))  - (VK)Sy(t). 

In  terms  of  the  new  time  variable  s defined  by  (16),  Theorem  1 reads  as  follows: 

Theorem  2 

Assume  that  either  a(t)€  L(0,+m)  or  h(() 6 L(0,  ■«■•).  Then  the  X force  will  be  annihilated 
in  finite  time  if  and  only  if 


y,  K \Z*cAs^- sris^y 


(27) 
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where  the  modified  time  variable  s is  given  by  (16),  and  Cx(s),  sx(s),  Cy(s),  and  sy(s)  denote  the 
normal-form  hyperbolic-like  GLF.  The  modified  parity-condition  parameter  Z*  is  unique  and 
given  by 


lim 


sx(s)  1 
cx(s)  Z»- 


(28) 


We  observe  that 


Q*  = KZ*, 


(29) 


and  liin,^{jv(s)/ci'(s)}  = Z*.  When  (27)  holds,  the  time  to  annihilate  X,  denoted  as  t/,  is 
determined  by  x(r«^)  = 0.  If  we  denote  the  quotient  of  the  two  normal-form  hyperbolic-like 
GLF  5x(x)  and  Cx(5)  as  i)x(f),  then  it  follows  from  (13)  that 


Vx(s(tJ‘))  = 


{xeCr(5o)  + y(iiy/{XK)IK)Sx{S(,)] 
{xoJrlso)  + yo(V(Aj»)/K)Cx(4o)}’ 


(30) 


where 

rtx(s)  = sx{s)lcx{s).  (31) 

5.  DETERMINATION  OF  THE  PARITY-CONDITION  PARAMETER 
We  will  now  show  how  knowledge  about  the  modified  parity-condition  parameter  Z*  for  one 
pair  of  attrition-rate  coefficients,  a(t)  and  bi(t),  allows  us  to  determine  Z*  for  a related  pair, 
a(()  and  bit)-  With  this  in  mind,  let  us  denote  Cx(5)  corresponding  to  a(t)  and  bU)  as 
Cx(5;  a,  b),  and  similarly  for  and  rix-  In  other  words,  we  will  now  write  (31)  corresponding  to 
the  attrition-rate  coefficients  a(t)  and  b(t)  as 


7jx(j;  a,  b)  = sx(s;  a,  b)/cx(s;  a,  b).  (32) 

In  this  notation,  we  will  write  (28)  as 

‘ lim  rix(s;a,b)^llZ*la,b].  (33) 


We  use  the  notation  Z*la,b]  to  show  that  the  modified  parity-condition  parameter  is  a 
functional  (i.e.  a function  for  which  the  independent  variables  themselves  are  functions),  which 
depends  on  only  the  attrition-rate  coefficients  a(t)  and  b(t).  In  other  words,  the  attrition-rate 
coefficients  are  functions  defined  for  t(,st<  + »,  and  the  parity-condition  parameter  depends 
on  these  entire  functions  (and  not  merely  particular  values  of  them). 

Our  main  result  is  Theorem  5,  which  gives  an  error  estimate  for  the  iq>proximation  that  we 
propose  for  Z*.  The  theoretical  basis  for  Theorem  5 is  given  by  Theorem  4,  which  (in  turn)  is  a 
consequence  of  Theorem  3.  The  proof  of  Theorem  3 follows  along  the  lines  of  well-known 
arguments  (see  p.  225  of  [26]). 

Theorem  3 

Comparison  Theorem:  Let  r(()  and  X|(()  satisfy 

with  initial  conditions 

x(t  - (•)  - a,  JCi(r  - 1#)  = a, 

{[l/a(/)J  dx/d/(/)},.^  - ft,  {(l/a(r)]  dx,/d/(0},.a  - ft, 

where  a(t)>0  and  bi(t)<b(t)foralll>tv.  Then  x,(0<x(l)foralll>ttas  long  as  x(t)>0. 
The  buk  theoretical  reauh  upon  which  our  numerical  determination  of  Z*  U based  is: 
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Theorem  4.  » i 

Assume  that  bt(t)  < b(t)  for  all  t > to-  Then 

Vx(s-,a,b)<\IZ*[a,b]<r,x(s;a,b)  + {(llZ*[a,bi])-vxU.a,o,)].  (34)  ^ 

Proof.  We  observe  that  [7]  r)xis',  o,  b)  satisfies  the  differential  equation  i 

dr,xlds(s  ■,a,b)  = \l{cx(s ; a,  b)]\  (35) 

with  ■nx{s-0;a,b)  = 0,  and  similarly  for  fix(.s,a,bi).  Theorem  3 (the  comparison  theorem) 
yields  that  Cx(s;  a,  b) > Cx(s\ a,  bi)  for  all  s > 0.  Thus,  for  all  s > 0 

drixlds(s;a,b)<drtxlis(s;a,bi), 

whence  integration  between  0 and  s yields  the  desired  result.  Q.E.D. 

Similar  to  the  observations  made  in  Remark  4 above,  we  observe  that  (33)  suggests  that  we 
estimate  Z*[a,  b]  with  Z defined  by 

i(#;a,h)=l/ijx(f;a.h),  (36) 


I 

i 

r 


where  s denotes  a suitably  chosen  value  for  s.  Moreover,  from  (35)  we  see  that  rix(s\  a,  b)  is  a 
strictly  increasing  function  of  s so  that  the  larger  we  take  s in  (36),  the  better  our  approximation 
becomes.  The  only  problem  (see  Remark  4)  is  that  a priori  we  do  not  know  how  large  to  take  i 
for  “satisfactory”  estimation  of  Z*.  Theorem  5,  however,  tells  us  exactly  how  large  to  take  I. 

Theorem  5 

Error  Estimate  for  Approximation:  Assume  that  bi(t)  < b(t)  for  all  t>t^  Let  /cff)  denote 
the  fractional  error  made  in  the  estimation  of  Z*[a,  b]  by  t(s  \a,  b),  i.e. 


a,b)-Z*[a,b] 

Z»[a,b] 


(37) 


tnen 

0<fB(i)<{m*la.  6,])-  ijx(f ; a,b,)}  ^(S-,  a,  b).  (38) 

Proof.  The  theorem  follows  by  simple  algebraic  manipulation  after  setting  s > f in  (34)  and 
using  (37).  Q.E.D. 

Thus,  we  have  presented  a method  for  numerically  determining  Z*ia,  h].  We  simply  pick  a 
large  value  for  s (we  denote  the  selected  value  as  I),  compute  SxiS)  and  Cx(S),  and  then 
compute  the  estimate  a,  b)  according  to  (36).  Theorem  5 allows  us  to  know  the  accuracy  of 
our  approximation,  which  can  be  improved  by  taking  f larger.  Thus,  we  can  numerically 
determine  Z*[a,  b]  to  any  specified  degree  of  accuracy  once  Z*[<i,  h|]  is  known.  In  the  next 
section  we  apply  this  theory  to  the  analysis  of  battles  modelled  with  offset  power  attrition-rate 
coefficients. 

6.  APPLICATION  OF  THEORY  TO  OFFSET  POWER  ATTRITION-RATE  COEFFICIENTS 

In  the  application  of  Theorems  4 and  5,  two  pairs  of  attrition-rate  coefficients  are  involved: 
one  pair  for  which  the  modified  parity-condition  parameter  is  known  [denoted  as  a(t)  and  hi(t)], 
and  one  for  which  it  is  to  be  determined  [denoted  as  a(t)  and  h(()].  Accordingly,  we  rewrite  (9) 
with  A>0as 

a(t)~k,(t  + Cy‘.  and  b(t)  * kt(t  + C + A)",  (39) 

where  (as  b^ore)  C2O.  We  will  refer  to  these  coefficients  (39)  for  which  A>0  as  potver 
attrition -rate  coefficients  with  “positive  of  set”.  If  we  choose 


► 


K~[uus  + iy^-'. 
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it  follows  from  (16)  that  the  modified  time  variable  s is  given  by 

siD^Min  + Df'U  + Cr*',  (41) 

and  the  invariant  J(s)  of  the  normal  form  (17)  simplifies  to 

/(j ; a,  h)  = /(j;  y,ii,p)  = (42) 

where  p- {ft + 1)11,  a = !/(/*  + 1),  p = (v  - ii)l{ii  + 1),  y = A • [Aj/(#t  + 1)]*®,  and  2 = /a  + 1/  + 2. 
Here  we  have  denoted  the  invariant  corresponding  to  the  attrition-rate  coefficients  a{t)  and  b{t) 
as  J{s ; y,  ft,  v),  since  we  may  take  y,  ft,  and  i/  as  a basis  for  generating  the  four  parameters  a,  p, 
y and  v that  explicitly  appear  in  the  right-hand  side  of  (42).  Furthermore,  we  will  denote  the 
normal-form  hyperbolic-like  GLF  that  correspond  to  J{s;y,tt,v)  as  Cx{s;y,fi,v)  and 
sxis',  y,  ft,  v). 

The  known  nsults[l]  that  we  use  in  the  Theorems  4 and  S are  for  the  case  of  power 
attrition-rate  coefficients  with  no  offset  [i.e.  set  A = 0 in  (9)] 

a(^)  = k,(^  + C)^  and  h,(t)  = Mr  + C)',  (43) 

where  C 2 0.  We  observe  that  f>i(r)  < b{t)  for  all  r > - C.  It  follows  that  J{s  ;a,bi)  = s^  and  [7] 

Z*[a.f»il  = P®’''’r(l-p)/r(p).  (44) 

Thus,  for  the  bound  on  Z*[a,  b]  = Z*(y,  ti,  v)  given  in  Theorem  4 and  the  error  estimate  for  our 
apT'Yoximation  (36)  given  in  Theorem  5,  we  have  [6] 

i,x(s:a.hi)  = p''-’'*T,(S),  (45) 

where  S = Ips'"^^,  q = l-p,  and  T„  denotes  a Lanchester-Clifford-Schlifli  (LCS)  function, 
which  is  analogous  to  the  hyperbolic  tangent  (see  Table  1).  These  functions  were  introduced  in 
[12]  and  redefined  for  reasons  of  force-annihilation  prediction  in  [6].  In  fact, 

t|(s ; a,  hi)  = tanh  s,  when  (i*v.  (46) 

This  result  is  one  of  our  reasons  for  introducing  the  normal  form  (7). 

We  have  thus  shown  that  the  following  theorem  holds. 

Tibk  1.  Lsnchwter-Cliffofd-SchlMi  functioat 

Kthtion  to  normat-fom  GLF 

fjt(*)  - F,(5)  txU)  - p"-*'*H,(S) 

Cr(*)-F,(S)  *r(*) 

when  l-p  aad 

StO-lpi'**" 

NOTE:  For  M • r,  m have 

(I)  exit)  - Cr(*)  ■ Fin(f ) ■ coih  t, 

(II)  »»(f ) - *r(*)  “ Hi/i(*)  - «Bh  I, 

(in)  a»(*)  “ nr{i)  ■ fintf ) - taah  i. 
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Theorem  6 

Assume  that  either  or  >/>  — I . Then  for  a battle  modelled  with  the  offset  power 

attrition-rate  coefficients  (39),  bounds  on  the  modified  parity-condition  parameter  Z*(y,  p,  u) 
are  given  for  y>0  by 

Vx(s-,  y,ii,  iix(s;  y,n,  + (47) 


where  q = l-p,  S = 2pi and  Tix(s;y,/i,v)  denotes  the  quotient  of  two  normal-form 
hyperbolic-like  GLF  for  the  attrition-rate  coefficients  (39),  i.e.  rix(s,y,p.,v)  = Sx(s‘,y,p„v)l 
cx(s\y,p.,v). 

It  follows  from  Theorem  6 (or,  equivalently.  Theorem  4)  that  if  we  approximate  Z*(y,  p.,  v) 
with  ^(s;  y,  p,  p)  defined  by 

Z{s;y,p,p)=\lrix(S-,y,p,p),  (48) 


then  bounds  on  the  fractional  error  made  in  this  estimate  are  given  by 


0</b(#)<P 


where  /sff)  denotes  the  fractional  error  and  is  defined  by  (37). 

The  right-hand  inequality  in  (49)  [equivalently,  (47)]  tells  us  exactly  how  large  to  take  i for 
the  estimation  of  Z*(y  > 0,  p,  v)  by  Z(f ; y,  p,  v)  to  any  specified  degree  of  accuracy.  The  LCS 
function  T,  is  involved  in  the  bound  on  the  fractional  error  /e(I)  in  this  estimate  when  p/  v. 
As  (46)  and  Table  1 show,  T,(5)  * tanh  s when  p = v.  Thus,  the  LCS  functions  as  redefined  by 
Taylor  and  Brown  [6]  yield  valuable  information  about  battles  modelled  with  not  only  the  power 
attrition-rate  coefficients  with  no  offset  (43)  but  also  the  offset  power  attrition-rate  coefficients 
(39).  Availability  of  tabulations  of  these  LCS  functions  is  discussed  in  [6]. 


7.  NUMERICAL  RESULTS 

In  this  section  we  will  examine  several  numerical  examples  to  show  how  the  modified 
parity-condition  parameter  Z*  may  be  numerically  determined  and  to  show  some  important 
insights  into  the  dynamics  of  combat  that  may  be  consequently  obtained.  In  order  to  numeric- 
ally determine  the  modified  parity-condition  parameter  for  the  offset  power  attrition-rate 
coefficients  (39),  we  must  use  knowledge  about  how  quickly  the  limiting  value  (i.e.  Z*[a,  bil)  of 
a hyperbolic-tangent-like  function  of  a related  pair  of  power  attrition-rate  coefficiente  with  “no 
offset”  (43)  is  reached  as  its  argument  is  increased  [recall  Theorem  6 and  (49)].  In  Fig.  3 we  see 


Hyptrbolic  - Lik* 
LCS  Function, 


238  3.  G.  Taylm  and  G.  G.  BkOWN 


Fif.  4.  Dependence  of  the  modified  parity-condition  parameter  Z*  on  the  modified  offset  parameter  y for 
the  offset  power  attrition-rate  coefficients.  The  modifM  offset  parameter  is  given  by  y = A ■ iXiUfi  1)]^, 
where  A is  the  offset  parameter  and£>/i-f'i''f2. 


that  this  limiting  value,  denoted  as  Z*{ji,  v)  = Z*[a,  hi],  is  quite  quickly  reached:  if  one  takes 
f = 10.0,  then  Z*(m.  v)  is  approximate  to  better  than  six  decimal  places  by  t(s,n,v)  = 
l/nx(s;  >').  where  ifx  is  given  by  (45).  Experimental  computing  for  various  values  of  m and  v 
and  comparison  with  the  known  value  (44)  for  Z*(m.  v)  bears  out  this  degree  of  accuracy  [i.e. 
speed  of  convergence  of  /(/;  /t,  v)  to  Z*]  for  essentially  all  allowable  values  of  n and  v.  Thus, 
Z(l;  /t,  v)  for  the  coefficients  (43)  has  essentially  converged  to  v)  when  s = 10.0,  and  by 
Theorem  6 or  (49)  we  know  that  the  same  is  true  for  i(S\  y,  n,  v)  for  the  coefficients  (39). 

We  have  accordingly  generated  by  this  procedure  the  results  shown  in  Fig.  4.  For  computing 
Vx  = ixicx,  we  have  used  the  series  solutions  shown  in  Tables  2 and  3.  (In  Tables  2 and  3 we 
have  for  convenience  denoted,  for  example,  sx(s ; y,  n,  v)  simply  as  sx(s ; M.  v),  i.e.  sx(s ; >') 
denotes  sx  corresponding  to  the  general  power  attrition-rate  coefficients  (9)  with  exponents  n 
and  V.)  The  series  were  obtained  by  solving  (17)  by  the  method  of  successive  approximations 
(see  [23]).  We  used  these  series  instead  of  developing  approximate  solutions  by  fiffite-difference 
methods  because  we  did  not  have  any  error  bounds  for  the  latter. 

Let  us  now  give  an  intuitive  interpretation  of  the  curves  shown  in  Fig.  4 of  the  modified 
parity-condition  parameter  Z*  plotted  vs  the  modified  offset  parameter  y.  In  Taylor  and 
Comstock  [7]  it  is  shown  that  Z*  may  be  considered  to  be  the  initial  Y force  level  that  leads  to  a 
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draw*  (i.e.  parity  between  the  forces)  in  the  battle  against  an  X force  of  “unit  strength" 


djt/dj«-y  with  x(j -0)*  1, 
dy/dj“-y(j)x  with  y(j  » 0)  = Z*, 


(50) 


where  /(s)  denotes  the  invariant  of  the  normal  form  (17).  Thus,  we  may  consider  Z*  to  be  “the 
Y equivalent  of  an  X force  of  unit  strength”  for  the  modified  battle  (50).  Now  let  us  consider 
the  general  power  attrition-rate  coefficients  (9)  with  exponents  ft  and  p.  As  we  did  in  Tables  2 
and  3,  we  will  denote  the  corresponding  / as  J(s\n,  p)  and  Z*(yt,  p)  to  stress  the  dependence 
on  /i  and  S'  (but  suppressing  that  on  y).  We  then  have  from  (42)  that  /(s;  1, 1)  = 1 -f  yly/(s)  and 
y(5;l,2)=- V(»Xl  + y/V(*))*.  From  (44)  we  find  that  Z*(l,l)*  1.000  and  Z*(l,2)-0.806  for 
y » 0.  Observing  that  fw  y > 1 we  have  /(* ; 1, 1) < /(s;  1. 2)  for  all  s 2 0,  it  is  intuitatively  clear 
from  (SO)  and  the  interpretation  of  Z*  as  a force  equivalent  that  we  must  have  Z*(l,  1)< 
Z*(l,  2)  for  all  y > 1 because  X always  has  greater  fire  effectiveness  against  Y when  y*  1 and 
p*2  than  when  /i » 1 and  v * 1.  However,  for  y near  zero,  the  situation  is  reversed  and 
Z*(l, 2)  must  lie  below  Z*(l,  1)  for  y near  zero.  Thus,  we  have  given  an  intuitive  explanation  of 
why  Z*(l,  2)  lies  below  Z*(l,  1)  for  y near  zero  but  above  it  for  y > 1 as  Pig.  4 shows. 

Next,  we  will  consider  numerical  results  for  a particular  battle  to  show  some  of  the 
important  insists  that  may  be  gained  into  the  dynamics  of  combat  from  our  new  results.  As  in 
(6, 12,23]  we  consider  S.  Bonder’s[3, 19]  model  (10)  for  the  constant-speed  attack  (ff  mobile 
forces  against  a static  defensive  position.  We  will  focus  on  the  new  results  given  in  this  paper 
(in  particuular,  the  prediction  of  battle  outcome  from  initial  conditions  without  explicitly 
computing  the  force-level  trajectories).  Input  data  and  computed  parameter  values  are  shown  in 
Table  4.  We  will  now  consider  two  cases;  G)  Jlo’*  1500 m;  and  (U)  Xp*  1250m. 

When  Xo*  1500m,  we  have  C«0  and  s«-0.  The  maximum  time  that  the  battle  can  last  is 
(mi  ~ 11.18  min,  since  at  this  time  the  advancing  attackers  overrun  the  defensive  position.  In 
this  case  Z*(y,  *')  ~ Z*(0.32, 1.1) « 1.381,  so  that  Theorem  2 teOs  us  that  JK  can  be  annihilated 

< > > > xb/y«  < 0.264.  By  (30)  the  X-force  annihilation  time  is  given  by  i|x(*(t«^))  ~ 2.739xo/yo. 
For  Xo  > 10  and  y«  > 50,  we  have  i)x(r«^)  * 0.54772  so  that  by  the  techniques  introduced  in  [6] 
we  find  ■ O.Vt.  These  computations  for  determining  s«^  involve  the  generation  of  a taM 
of  tx,  Cx,  and  ijx  f«r  y“0.32,  ^">'*1  (*««  (8D-  Hence,  (36)  yields  t/«  10.25  min  and 
r/  « 125.7  m.  Flvther  results  are  given  in  Table  5. 

When  Xp  » 1250  m (see  Fig.  S of  [12]),  we  have  C » 1.864  min,  sp = 0.0255  and  * 9.32  min. 


*Ia  other  trotdi,  s(i)  ted  r<f ) > 0 for  d « e (0,  'f  ■)  hut  Mb - 0 - fim .^y(t ). 
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Table  4.  Particulars  for  the  numerical  examples 

Table  5.  Annihilation  of  Uie  X force  as  a function 

1.  Input  data 

= •>=  1 

(Xa/y«) 

t/(min) 

r/(m) 

a,  = 0.06X  casualties/min/y  unit 

AixO.by  casualties/min/X  unit 

0.250 

14.09 

K.  = 1500  m.  = 2000  m 

0.200 

10.25 

125.7 

V • 5 miles/hr 

0.167 

8.80 

319.4 

2.  Panmtier  values 

1 1.18  min  and  x/  = x(r  = 0) 

= 2.48. 

k.  = 5.364  X I0~’X  casualties/min/V'  unit 
tt  = 4.023  X I0~’y  casualtie$/min/3f  unit 

p = <1  = 112 

A = 3.728  min.  y = 0.320  (casualties/min)''^ 

Tabic  6.  Annihilation  of  the  X force  as  a function 
of  the  initial  force  ratio  for  Jlo  - 1250  m 


(Xo/yo) 

l,*(min) 

r/(m) 

0.250 

10.87 

J 

0.200 

8.17 

154.4 

0.167 

6.93 

320.4 

= 9.32  min  and  .19  = x(r  * 0)  = 1.74. 


In  this  case  X can  be  annihilated  < = = > xo/yo<  0.281  with  the  X-force  annihilation  time  given 
by  = (1.001«o+0.009)/(0.127iio+0.366),  where  Uo=  Wyo-  Numerical  results  are  given  in 

Table  6.  Finally,  these  parametric  results  should  be  contrasted  to  those  previously  possible  (e.g. 
compare  them  with,  for  example,  the  sinjgle  force-level  trajectOTy  for  R0  = 2000  m shown  in  Fig. 
5 of  [12J. 


8.  DISCUSSION 

S.  Bonder  (3, 18, 19]  has  emphasized  the  shortcomings  of  constant-coefficient  Lanchester- 
type  combat  models.  Work  by  Bonder [3, 18],  Clark[2S],  and  others[4]  on  the  prediction  of 
Lanchesler  attrition-rate  coefficients  (see  Taylor  and  Brown  [12]  for  further  discussion  and 
references)  has  generated  interest  in  variable-coefficient  models.  Moreover,  there  is  not  only 
intrinsic  interest  (see  [3. 19])  in  the  model  (7)  but  also  interest  for  obtaining  insights  into  the 
behavior  of  complex  Lanchester-type  system  models  (for  example,  the  Bonder-IUA  model  (see 
[4,3,14]))  that  have  been  enriched  in  military  detail  see  [4,3,14-17]).  The  attrition-rate 
coefficients  in  (7)  represent  the  fire  effectiveness  of  the  combatants  and  allow  us  to  model 
temporal  variations  in  fire  effectiveness  on  the  battlefield.  Interest  in  the  general  power 
attrition-rate  coefficients  (9)  is  provided  by  S.  Bonder’s  [3, 3, 19]  constant-speed  attack  model* 
(10),  (II)  and  his  examination  the  range  dependence  of  attrition-rate  coefficients  for  various 
weapon  systems  (tee  pp.  196-200  of  [4]). 

We  have  given  resukt  that  allow  one  to  study  the  variable-coefficient  model  (7)  (especially 
with  the  general  power  attrition-rate  coefficienu  (9))  afanott  as  easily  and  thorou^y  as 
Lanchester’t  clastic  constant-coefficient  model  (1).  Taylor  and  Comstock[7]  (see  Theorems  I 
and  2 above)  have  shown  how  to  predict  force  annihitotion  without  having  to  spend  the  time 
and  effort  of  explicitly  computing  force-level  trajectories.  Using  their  theoretical  results,  we 
gave  results  in  a previous  paper[6]  that  made  combat  modelled  by  power  attrition-rate 
coefficienu  with  no  offset]  (Le.  i4>«0  in  (9)]  almost  as  easy  to  analyze  u the  constant- 
coefficient  case.  The  resuhs  ^ the  paper  at  hand  allow  one  to  analyze  combat  modelled  by 
power  attrition-rate  coefficienu  with  positive  offsets  [Le.  A >0  in  (9)]  just  as  conveniently. 

Theorem  1 (see  also  Theorem  2)  is  the  generalization  of  the  clastic  constant-coefficient 
result  (3)  to  cates  of  time-dependent  attrition-rate  coefficienu.  However,  one  needs  to  know  the 

the  ranee  between  flrer  and  tnrtet  chnneee  dwins  the  emmement. 

tModellins.  for  eunvle,  combat  between  two  weapon  lyatema  with  the  lanw  maximuffl  effective  ranae. 

tModeUns.  for  exam^,  combat  between  two  weapon  iyateme  with  dWerant  maximuffl  effective  ransea. 
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value  of  the  so-called  parity-condition  parameter  Q*  in  order  to  predict  force  annihilation  in 
specific  instances.  In  this  paper  we  have  presented  theoretical  considerations  (see  Section  5 
above)  for  the  noniterative  numerical  determination  of  the  parity-condition  parameter.  We 
applied  our  general  theory  to  the  specific  case  of  general  power  attrition-rate  coefficients  (9) 
(see  Section  6)  and  illustrated  these  theoretical  results  by  considering  some  numerical  examples 
(see  Section  7). 

Curves  of  the  modified  parity-condition  parameter  Z*  plotted  against  the  modified  offset 
parameter  y such  as  those  shown  in  Fig.  4 allow  one  to  analyze  parametricaUy  “modem" 
combat  modelled  with  the  general  power  attrition-rate  coefficients  (9).  For  example,  we  can  now 
parametrically  (for  example,  varying  the  maximum  effective  range  of  the  defender’s  weapons) 
determine  whether  the  defender  will  be  overrun  in  Bonder’s(3,S,  19]  constant-speed-attack 
model  (10)  with  attrition-rate  coefficients  (II)  without  having  to  compute  the  entire  force-level 
trajectories.  We  have  illustrated  this  analysis  capability  with  some  numerical  examples,  which 
show  that  the  defender’s  annihilation  (i.e.  saturation  of  his  defensive  position  with  offensive 
fire)  depends  on  the  initial  force  ratio  (of  defender  to  attacker)  being  below  a certain  threshold 
value.  Our  new  results  allow  one  n6U>nly  to  determine  easily  such  force-ratio  thresholds  of 
survivability  but  also  to  study  their  dependence  on  weapon-system-capability  parameters. 

Our  new  results  let  us  conveniently  obtain  much^ valuable  information  about  the  model  (7). 
The  classic  ordinary  differential  equation  theories  (see,  for  example,  Ince[26]  were  inadequate 
to  answer  many  important  questions  (for  example,  “who  will  win?  Be  annihilated?’’)  about  such 
combat  models.  Previously  one  was  limited  to  only  computing  force-level  trajectories,  but  now 
we  can  predict  battle  outcome  (in  particular,  force  annihilation)  without  explicitly  computing 
force-level  trajectories.  Moreover,  these  new  results  facilitate  parametric  analysis  of  such 
combat  situations.  S.  Bonder  [22]  has  suggested  that  an  increased  emphasis  be  place  on 
parametric  analyses  in  systems  analysis  studies  (see  pp.  21-22  of  [22]).  In  particular.  Theorems 
1 and  2 explicitly  exhibit  a tradeoff  between  quality  (as  quantified  by  the  relative-fire- 
effectiveness  parameter  A|>  and  the  parity-condition  parameter  (?*)  and  quantity  (as  quantified 
by  the  initial  force  ratio  Xo/yo)  of  two  weapon  systems  in  combat  against  each  other.  In  other 
words,  one  can  use  an  expression  like  (21)  to  develop  quantitative  insights  into  how  the  quality 
of  a weapon  system  may  be  substituted  for  sheer  numbers.  Moreover,  an  unanswered 
theoretical  question  is  to  determine  how  the  parity-condition  parameter  Q*  depends  on  the 
combat-intensity  parameter  A/  and  the  relative-fire-effectiveness  parameter  A*.  Finally,  our 
results  here  are  signposts  as  to  the  difficulty  of  analytically  extracting  information  (particularly 
parametric  information  without  excessive  computations)  from  variable-coefficient  I,anchester- 
type  models  such  as  (7). 
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